Precise toppling balance, quenched disorder, and universality for sandpiles 
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A single sandpile model with quenched random toppling matrices captures the crucial features of 
different models of self-organized criticality. With symmetric matrices avalanche statistics falls in 
the multiscaling BTW universality class. In the asymmetric case the simple scaling of the Manna 
model is observed. The presence or absence of a precise toppling balance between the amount of 
sand released by a toppling site and the total quantity the same site receives when all its neighbors 
topple once determines the appropriate universality class. 
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Self-organized criticality (SOC) occurs in the non- 
linear transport of some physical entity, like energy, sand, 
or stress, through a system of linear size L under a con- 
stant, slow external drive Q. 0.0.0. The transport has 
breakdown features, with intermittent bursts of activity 
called avalanches. At long times the probability distribu- 
tion of avalanche sizes becomes critical since they do not 
reveal characteristic scales between L and the minimal 
length specified in the model. In spite of many efforts, 
the physical mechanisms underlying the different forms 
of scaling realized in such models and the related univer- 
sality issues remain poorly understood. The probability 
distribution of the number of topplings in an avalanche of 
the Bak, Tang and Wiesenfeld (BTW) @ sandpile, his- 
torically the prototype of SOC, has been shown recently 
to violate the finite size scaling (FSS) ansatz, assumed 
for many years, and to obey a peculiar form of multiscal- 
ing 0,0- On the other hand, in the Manna stochastic 
sandpile |al the corresponding distribution is known to 
obey FSSlaGiJIIIl- 

Understanding the key mechanisms 
at the basis of these radically different forms of scaling 
remains a major challenge which should be faced also in 
the perspective of new concrete applications . 

In this Letter we study sandpile models 0] on lattices 
with quenched disorder. This means that different bonds 
of the lattice allow flow of different, but constant num- 
bers of sand grains through them. We further distinguish 
between an undirected case, in which the flow is identi- 
cal in the two directions for each bond, and a directed 
one, in which this flow can be asymmetric. We show 
that, once averaged over the possible realizations of the 
quenched disorder, the avalanche size distribution of the 
model with symmetric flow in each bond obeys the same 
multiscaling as the BTW model, while asymmetry in the 
flow leads to FSS with the exponents of the Manna sand- 
pile. Thus, for a disordered sandpile the two main SOC 
universality classes can both be realized simply upon en- 
forcing or releasing a local symmetry requirement. 

The deterministic sandpile model on a square lattice 
of size L is described usingan integer 'toppling matrix' 
(TM) A of size L 2 x L 2 |. At lattice site k there is 
a column of sand grains. The system is externally 



driven by adding a single grain at a time at a randomly 
selected site i: hi — > hi + 1. Each site has a threshold 
height Hi of stability. If hi > Hi, the sand column at i 
topples and grains are distributed to other sites. Conse- 
quently, all heights are updated as: hj — * hj — Ay where 
An > and Ay < for i ^ j. The threshold heights 
are chosen as Hi = An. Grain conservation is assured 
by putting Ajj = — ^<y ■ Grains must flow out of 
the system through the boundary sites to maintain a sta- 
tionary state. The BTW model is a special case of this 
formulation where A^ = 4, Ay = — 1 for \i — j\ = 1, and 
Aij = for \i — j\ > 1 |5(. In contrast, in the stochas- 
tic sandpile model || , annealed randomness enters in the 
grain distribution process upon toppling. Indeed, in this 
case each grain of the toppling site is transferred to a 
randomly selected neighboring site. For both BTW and 
stochastic sandpile grain number is conserved, bound- 
aries are open, and, when a site topples grains are dis- 
tributed to its neighborhood. In spite of these common 
features and similarities [ljj, the difference in behavior 
of the two models concerns even the type of scaling. 

In the Manna sandpile, the randomness in the choice 
of the neighbors getting a grain from the toppling site 
can be regarded as an "annealed" disorder. Indeed, the 
random choice of the TM elements is constantly updated 
during the course of a given avalanche. In contrast, a 
"quenched" disorder is realized in the models discussed 
here, for which a random realization of A is maintained 
for full samples of avalanches, and an average over the 
distributions of many samples created is performed even- 
tually. 

Our A^ is nonzero only for i — j or \i — j\ = 1. 
However, unlike in the BTW model, Ay takes random 
values. These are chosen among the negative integers 
— 1, —2, .., —to when \i — j\ — 1. Thus, when one of the 
sites joined by a nearest neighbor bond topples, this bond 
can allow flow of more than one grain, unlike in the BTW 
case. When the TM is symmetric we call the model undi- 
rected and assign only one random integer Ay = Aji in- 
dependently to each nearest neighbor bond (Fig. 1(a)). 
In the directed case the TM is asymmetric and both Ay 
and Aji are assigned independently drawn random inte- 
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FIG. 1: Random grain flow distribution along the bonds of 
a 3 x 3 square lattice: (a) undirected case, (b) directed case. 

gers (Fig. 1(b)). In this way a toppling at one end of a 
bond may send, through the same bond, a different num- 
ber of grains from that at the other end. At any site i 
we put Hi = An = — Yjj-^t A-ij ■ The total number of 
grains received by the site i when all neighboring sites 
topple once is H[ — — Sj^jAjj. We call Hi and H[ the 
out-degree and in-degree of site i, respectively. 

Irrespective of whether A is symmetric or not, the 
models enjoy the Abelian property of the BTW sand- 
pile 0| . This allows to establish a number of general 
exact results concerning the stationary state, recurrent 
configurations, etc., which hold in both the directed and 
the undirected case Q- An important notion is that 
an avalanche can be decomposed into waves of toppling 
[13. Ilfij . To this purpose one considers the sequence of 
topplings following the addition of the seed grain at site 
0. The first wave is the set of all topplings that follows 
the first toppling at while itself is prevented from a 
possible second toppling. If is still unstable after the 
first wave, the second wave starts by allowing a second 
toppling at 0, and so on. 

A peculiar property of waves with a symmetric TM is 
that the set of lattice sites which topple has no holes. 
For example, a single untoppled site can not be fully 
surrounded by toppled sites. Indeed, in the undirected 
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FIG. 2: Multiply toppled sites within avalanches are shown 
by circles of different colors: 1 (black), 2 (red), 3 (blue), 4 
(green), 5 (magenta), 6 (orange) for the (a) undirected model 
(b) directed model. 

model the equality Hi = H[ is strictly maintained at all 
sites except at the boundary, which implies that a site 
must topple irrespective of its height, if all its neighbors 
topple once. Thus, waves in the undirected model, like 
in the BTW model, have no holes. In addition one can 
show that all sites involved by a wave topple only once, 
like the seed site. All this is not true anymore with asym- 
metric TM. In this case the waves can have holes and 
can include sites which topple more than once. Indeed, 
in the directed model, H t ^ H[ in general and a site with 
Hi > H[ does not topple even if all its neighbors topple 
once. This creates a single site hole in the avalanche. 
On the other hand, if a site i has out-degree sufficiently 
smaller than its in-degree, it may topple for the second 
time even if none of its neighbors have toppled for the 
second time. Thus, in general different sites belonging to 
a given wave topple a different number of times. 

The compactness and uniformity of waves in the undi- 
rected model leads us to expect for them and for the 
resulting avalanches a structure similar to that of the 
BTW model (Fig. 2(a)). The non-uniformity introduced 
by disorder should not be relevant at large scales when 
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FIG. 3: (a) Plot of <r{q) vs. q for the BTW (solid), undi- 
rected (dotted), Manna (dashed), directed (dot-dashed) mod- 
els. Comparison of da(q)/dq vs. q plots between (b) BTW 
(solid lines) and undirected (dotted lines) models and for the 
(c) Manna (solid lines) and directed (dotted lines) models. 



one counts toppling events. The situation is quite differ- 
ent for the directed model: the structures of waves and 
of avalanches in this case (Fig. 2( b)) look in fact similar 
to those of the Manna model where the numbers of 
grains transmitted upon in the two directions of a given 
bond are different in general, and this imbalance is main- 
tained dynamically due to the stochastic distribution of 
sand grains. 

In the BTW model, the probability distribution 
Prob(s, L) of the total number of topplings, s, in an 
avalanche has been found recently to obey a multiscaling 
ansatz 7]. On the other hand, it is pretty well estab- 
lished [10( by now that in the Manna stochastic sandpile 
this distribution obeys simple FSS as: 



Prob(s, L) ~ s 



(1) 



where the scaling function f(x) ~ constant in the limit 
of x — > and f(x) approaches zero very fast for x » 1. 
The exponent r and the dimension D fully character- 
ize the scaling of Prob in this case. One immediate 
way to check validity of Eqn. (1) is to attempt a data 
collapse by plotting s r Prob vs. s/L D with trial values 
of the exponents. We collected extensive data for both 
our models [m = 4) for L = 128, 256, 512, 1024, 2048 
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FIG. 4: Autocorrelation function of the wave time series of 
the undirected and the directed (inset) sandpile for L=128, 
256 and 512. 



and 4096, namely 50 million avalanches in 500 indepen- 
dent configurations for L = 128 down to w 1.1 million 
avalanches for 9 configurations for L = 4096. The first 
~ 4L 2 avalanches were skipped to reach the steady state. 
For the directed case collapse works very nicely giving 
r w 1.28 and D w 2.75, close to the most reliable es- 
timates of the Manna sandpile exponents ^j. For the 
undirected model the collapse does not work for a single 
set of t and D and for all values of s and L. This is found 
also in the BTW sandpile. 

A more reliable and quantitative check of the valid- 
ity, or violation, of FSS, is based on the evaluation of 
the various moments of Prob 0, 0> Ell The q-th mo- 
ment is defined as (s q ) = j s 9 Prob(s, L)ds. Assuming 
that FSS holds, it is easy to show that (s q ) ~ L"^ with 
the moment exponent given by <j(q) = D(q — r + 1) for 
q > t — 1 and a(q) = for < q < r — 1. In the case 
of multiscaling a should have a nonlinear dependence on 
q. A comparison of o~(q) for two given models is also a 
key to establish if they belong to the same universality 
class, or not. The value of cr(g) was determined from 
the slope of the plot of \og(s q (L)) vs. logL for L = 
1024, 2048 and 4096 with an error w 0.01 and for 251 
values of q between and 5 (Fig. 3(a)). The derivative 
of a is determined by the finite difference method. Slow 
but monotonic increase of da(q)/dq with q clearly indi- 
cates the multiscaling in BTW as well as in undirected 
models (Fig. 3(b)) whereas a saturation of da(q)/dq in- 
dicates validity of FSS for Manna and in directed models 
(Fig. 3(c)). To measure deviations quantitatively we 

dpfmp a nnantitv X , — 9 L 'dq) a -(da(q) / >dq) b \ t, • 

aenne a quantity A a , b - z {dcr{q) / dq)a+ ( da ( q )/ dq)b ■ « is 
observed that after some initial fluctuations XBTW,undir 
has a maximum value 0.33% within q=2 and 5 whereas 
XBTW,direc gradually increases to 5.5 % at q =5. This 
analysis implies that the undirected model is almost neg- 
ligibly different from the BTW model in the da(q)/dq 
vs. q plot whereas the deviation of directed model from 
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the BTW model is much larger and gradually increases 
with q. Similarly XManna,direc is limited witin 0.91 % 
whereas XM a nna,undir gradually increases to 6.5 % at q 
=5 which also implies that the directed model behaves 
very similarly to the Manna model, and is very much dif- 
ferent from the BTW model. The above described results 
concerning X for the various couples of models are alto- 
gether strongly supporting the conclusion that while the 
directed model belongs most likely to the Manna uni- 
versality class, the undirected one has the multiscaling 
features known to be peculiar of the standard BTW sand- 
pile. 

By decomposing a large sequence of successive 
avalanches into waves in the undirected and directed 
cases, we obtained global wave size distributions which 
obey FSS with the exponents expected for the BTW 
model and the Manna model ^||, respectively. For 
the BTW sandpile globally sampled waves have a size 
distribution with a form as in Eq. (1), with r w = 1 and 
D w = 2. In the case of the Manna stochastic sandpile 
waves can not be defined as in the deterministic Abelian 
sandpiles, but a wave-like decomposition was proposed 
in Ref. |l6(. The global size distribution of the corre- 
sponding waves obeys FSS with the same exponents ob- 
tained for the avalanche distribution |l6j | . As already re- 
marked above waves can be consistently defined 0] in the 
same way for each quenched disorder realization of our 
directed and undirected models. The global wave scal- 
ings obtained here further support the expectation that 
they fall in the Manna and BTW universality classes, 
respectively. 

Further insight into the different behaviors of the di- 
rected and undirected models can be obtained by ana- 
lyzing the wave time series j_si , 82,83,..} of the sizes of 
successive waves as in ref. 16]. In Fig. 4 we plot the 
autocorrelation function: 
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where the expectation values refer to samples with dif- 
ferent L and include quenched disorder averaging. The 
plots in Fig. 4 are fully consistent with similar ones 
for the BTW and Manna sandpiles 16]. While in the 
directed case the autocorrelation function is essentially 
zero as soon as t > 0, in the undirected model it grows 
steadily with L, and approximately scales as C(t,L) ~ 
t- T 'g(t/L D ') with t c w 0.35 and D c w 1. These expo- 
nents should be compared to 0.40 and 1.02, respectively, 
as determined for the BTW model 0]. This long range 
autocorrelation must be a consequence of the coherent 
and uniform spatial structure of each wave in the undi- 
rected case. In the directed model correlations are de- 
stroyed by the much more irregular pattern of topplings, 
with inhomogeneities and holes, in each wave. The cor- 



relation patterns show marked self-averaging, being re- 
producible on the basis of very few disorder realizations. 

The local out/in degree balance Hi = H[ at all sites 
in the undirected model is essential for the BTW multi- 
scaling behavior to prevail. Numerically, with quenched 
disorder realized as described above, we find that the 
density of unbalanced sites with Hi ^ H[ in the directed 
model is around 0.88 and those of the sites with H[ > Hi 
and Hi < Hi are equal to 0.44. Now we ask if there 
is any critical density of unbalanced sites which demar- 
cates the behaviors of the undirected and directed mod- 
els. To study this we first generated an asymmetric TM 
in which the fraction of the bonds with unequal A values 
is found to be « 0.75. We tuned this fraction, and thus 
the density of unbalanced sites, by randomly selecting 
these bonds and making their A values equal by assign- 
ing a random integer number between —1 and —to. We 
find that even the presence of as low as 5 percent bonds 
with unequal A values is sufficient to destroy the mutis- 
caling and to ensure FSS as in the directed model. Thus, 
as soon as precise toppling balance is broken, FSS holds, 
and the universality class turns into that of the Manna 
sandpile. The transition to Manna behavior does not 
require a nonzero threshold density of unbalanced sites. 
Thus, for the undirected model the symmetry of precise 
toppling balance is a crucial requisite for the multiscaling 
to hold. This requisite is of course satisfied also by the 
ordinary BTW sandpile. 

To conclude, we studied sandpile models with 
quenched disorder where the elements of the TM are ran- 
domly assigned. With asymmetric TM the precise top- 
pling balance between in- and out-degrees at each site 
is not maintained. This imbalance suppresses the wave 
correlations leading to the BTW-like multiscaling behav- 
ior of the avalanche size distribution and results a FSS 
regime in the universality class of the Manna stochas- 
tic sandpile. Thus, a symmetry mechanism underlies the 
puzzling difference between BTW and Manna scalings. 
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